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It is shown, by simulation and theory, that circularly or elliptically polarized terahertz radiation 
can be generated when a static magnetic (B) field is imposed on a gas target along the propagation 
direction of a two-color laser driver. The radiation frequency is determined by yCf-l-tuI/d-|-a;c/2, 
where ujp is the plasma frequency and cjc is the electron cyclotron frequency. With the increase of 
the B field, the radiation changes from a single-cycle broadband waveform to a continuous narrow- 
band emission. In high-B-field cases, the radiation strength is proportional to tUp/wc. The B field 
provides a tunability in the radiation frequency, spectrum width, and field strength. 


Terahertz (THz) spectroscopy and coherent control 
have been wid^ applied in physics dHil, biology @ 


and medicine [6|. These applications can potentially 
benefit from THz radiation sources from gas (rl-fl^ or 
solid 2^, 211 plasmas irradiated by fs intense laser pulses 
thanks to their high radiation strength and bandwidth 
up to 100 THz. Recently, powerful THz radiation of 
multi-MV/cm [^, has been efficiently generated via 

a two-color laser scheme in which a fundamental pump 
laser is mixed with its second harmonic in gases Q. Basi¬ 
cally, such radiation generated by linearly-polarized laser 
drivers is linearly polarized although, in some conditions, 
the linear polarization becomes elliptical during propa¬ 
gation due to modulation of the laser phase and polar¬ 
ization in gas plasma (^ . To achieve radiation with 
controllable polarization to further broaden the THz ap- 



In this Letter, we propose a scheme in which a static 
B field is imposed along the propagation direction of 
a two-color linearly polarized laser driver to generate 
narrow-band THz radiation of circular or elliptical po¬ 
larization with the relative phase between the two radi¬ 
ation field components fixed at 7r/2. The radiation rota¬ 
tion direction can be controlled by the B-field sign. At a 
field strength of lOOT, the electron cyclotron frequency 
Wc = eBolmeC is much higher than the plasma oscil¬ 
lation frequency ujp = (rif. is the formed 

plasma density), so the radiation frequency is almost at 
Wc, and therefore, it can be smoothly tuned by the B- 
field strength Bq. In this case, the B field dominates 
over the plasma oscillation and the former governs the 
trajectory of the plasma electrons, which causes such ra¬ 


diation properties. Because of Wc ujp, the radiation has 
a many-cycle waveform rather than a single-cycle wave¬ 


form [7|, I9l-ll2l. l28l| in the case without the B field. Thus, 


the current radiation has a narrow-band spectrum. 


Magnetic fields at tens of teslas are widely available 
in the form of dc or ms-pulsed, nondestructive mag¬ 


nets [3l|, where the highest one reaches lOOT. Via de¬ 
structive methods, 600T ^s-pulsed B-fields were avail¬ 
able more than a decade ago [H, 331. Nanosecond- 


laser-driven capacitor-coil experiments demonstrated ns- 
pulsed B fields of 1500T recently 3^, which also has 
significant applications in novel magnetically assisted in¬ 
ertial confinement fusion 


We first demonstrate the scheme sketched above 
through particle-in-cell (PIC) simulations with the two- 
dimensional (2D) KLAPS code [s^, in which the field 
ionization of gases is included. The pump laser wave¬ 
length is fixed at l^m (or the period tq = 3.33 fs) and the 
second laser frequency is at the second harmonic of the 
pump one. The two pulses propagate along the -|-x direc¬ 
tion with linear polarization along the z direction. They 
have the same spot radius rg = 150/rm and duration 50fs 
at full width at half maximum. Peak intensity of the 
pump pulse is 2 X 10^® W/cm^ with the energy 42mJ and 
the second pulse has the peak intensity 5 x 10^^ W/cm^ 
and energy llmJ. A helium gas slab is taken with a uni¬ 
form density 1.22 x 10^® cm“^ (the corresponding plasma 
frequency ujp = ITHz after the complete first-order ion¬ 
ization by the used laser pulses [l^) and a length 320/rm. 
The resolutions along the x and y directions are 0.01/rm 
and 0.25/rm, respectively. Initially, in the gas region, 
four simulation particles per cell denoting gas atoms are 
adopted. 


Figure 1 shows spatial distributions of the THz radi¬ 
ation propagating along the -x direction in the vacuum. 















2 



X (nm) * (H-*®) 


FIG. 1. [(a) and (b)] Snapshots of the THz electric fields 

(MV/cm) at the time of 0.7 ps and (c) the field distributions 
on the axis {y = 0) at 2 ps, where an external B field of 178T 
is imposed. The broken lines in (c) correspond to the case 
without the B field. 


which is generated with an external static B field of 178T 
imposed along the +x direction. As a comparison, the 
radiation generated without the B field is also displayed 
by the broken lines in Fig. 1(c), illustrating that the ra¬ 
diation has only the z-direction component and a near 
single-cycle waveform, as shown in previous experiments 
and simulations 0) ffllll- With the B field, the radia¬ 
tion also has the y-direction component in addition to 
the z-direction one. The two components have the same 
frequency, higher than that in the case without the B 
field, and a constant phase displacement. 

To further analyze the radiation properties, we take 
the temporal waveform observed in the left vacuum 5 ^m 
away from the vacuum-gas boundary, as illustrated in 
Fig. 2. The low-frequency part below 1 THz has been 
filtered and the high-frequency part is retained. One 
can see the radiation frequency at 5 THz, equal to the 
cyclotron frequency ujc- This frequency deviates from 
the central frequency around the plasma frequency Wp = 
ITHz of the radiation without the B field. 

Without the low-frequency or dc part, the two compo¬ 
nents Ey and show nearly the same strength and a 
constant phase displacement of 7r/2, i.e., circular polar¬ 
ization. When we reverse the B field to -x direction, the 
phase displacement is changed to —7r/2, i.e., the rotation 
of the CP radiation is also reversed, as observed in Fig 
2. Note that, in real applications, the dc part of the ra¬ 
diation could be rapidly absorbed as soon as it touches a 
material rather than the vacuum. The effective radiation 
interacting with a sample should be a CP wave as shown 
in Fig. 2(a). 



FIG. 2. (a) Temporal waveforms of the THz electric fields on 
the axis and (b) the corresponding spectra, where an external 
B field of 178T is imposed. The components below 1 THz are 
filtered. 


Now, we explain the radiation observed. The THz ra¬ 
diation process [Ulllll takes place as follows: first, a net 
current and plasma are formed via ionization, the current 
drives an electrostatic oscillating field in the plasma, and 
then this field is converted into electromagnetic radiation 
at the plasma boundaries. Without an external B field, 
electrons released from atoms have velocities only along 
the laser polarization, say the z direction, and therefore 
the generated radiation is also linearly polarized along 
the z direction. With the B field along the x direction, 
the electrons rotate in the y-z plane and then have veloc¬ 
ities in both the y and z directions. Hence, the radiation 
has the components along both the y and z directions. 

Frequency and waveform.—We set Ey and E^ as the ra¬ 
diation or oscillation electric fields formed in plasma. The 
nonrelativistic motion equation of an electron is dvy/dt = 
—cEylnif. — uicVz and dv^/dt = —eEzIme + iVcVy, where 
both the external B field and the radiation fields much 
below relativistic strengthes are considered. One easily 
derives 


dV 

dt 


eS . 

- lOJcV, 

me 


( 1 ) 


where V = Vy + iv^, £ = Ey + iE^, and i is the imaginary 
unit. According to the wave equation (V^—9^/c^i9f^)E = 
— (A'Kenel(?)d'vldt in plasma with a density Ue, one ob¬ 
tains 



dTrerie dV 
dt 


( 2 ) 


We first use Eq. (I2|) omitting the spatial differential term 
to look for the oscillation frequencies of the radiation 
source at a given position. Set V and £ with a frequency 
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oj. According to Eqs. o and m, one obtains: 


w± 



( 3 ) 


Note that one can derive uj = ujp from Eqs. m and m in 
the same way for Bq = 0. The three frequencies Wp and 
uj± correspond to the cutoff frequencies for wave prop¬ 
agation in unmagnetized and magnetized cold plasmas, 


respectively 


The B field separates the oscillation 


frequency from ujp into two frequencies: one above ojp 
and the other below it. In particular, when ojc S> Wp, ui+ 
approaches Wc, which provides a robust method to con¬ 
trol the radiation frequency by the B-field strength. The 
simulation results in Eig. 2 and the red line in Eig. 3(b) 
are in good agreement with Eq. The theoretical val¬ 
ues are 5.19 THz and 0.19 THz, compared to 5 THz and 
0 THz in the simulations with the numerical resolution 
of 0.5 THz. 

Performing Eourier transform to Eqs. o and m, 
one obtains the dispersion relation of the radiation wave 
along the -x direction, which has the refractive index 
77= — (yP- -I- wwc). Under the condition of uj+ ~ 


ujc 3> ujp, ry ~ 1, and therefore, the a;+ component of 
the radiation generated even in deep plasma can propa¬ 
gate to the vacuum with little attenuation [s^. Hence, 
the radiation is many-cycle and narrow-band as shown 
in Figs. 1 and 2. This is different from the single-cycle 
and broadband radiation without the B field because of 
its central frequency at Wp and ry = 0. 

Polarization.—Since high-B fields are required by 
frequency-tunable radiation, we consider such B fields 
as 


UJl^ UJci^ ujp, 


( 4 ) 


where ujl is the laser fundamental frequency. With 
Wc 3> Wp, the B field dominates over the plasma os¬ 
cillation, and thus, the velocity of an electron satisfies 
Vyj = Vj cos{uJct -I- Oj) and Vzj = Vj sin(a;ct + Oj). With 
wl ^ Wc the initial phase 9j can be considered as roughly 
the same for all electrons, since they are released only at 
the laser peak within a few cycles 18| when 50fs laser 
duration is used here. Then the average electron velocity 
just after the passage of the pulses can be written by 


f VyO = VQCOs{uJc.t + 9o), 
\ Vzo = no sin(a;ct -b 0o). 


( 5 ) 


We replace the electric fields with the vector potentials 
Ay and Az in Eqs. m and ©■ From the two equa¬ 
tions, one obtains dV/dt = {e/mec)dA/dt — iuJcV and 
—d"^/c^dP)A = {ATTene/c)V, respectively, where A = 
Ay + iAz . Here we are interested in the higher-frequency 
component with uj+ Wp. This component with ry cs 1 
can propagate in the plasma as in a vacuum, which could 
be considered as a plane wave. Therefore, the corre¬ 
sponding electron velocity follows —iw+V ~ dV/dt and 


—fWcV ~ {uJc/uj+)dV/dt, with which the motion equa¬ 
tion is rewritten by (1 — a;c/a;+)(V — Vo) = (e/mec)A. 
Inserting this equation of motion into the wave equation 
expressed by A, one obtains 


- 
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A dP 0^(1 — ujc/uj+) 


4Trene , . 
A = - Vo, (6) 


where Vo is given by Eq. ®. Equation ([5]) describes ra¬ 
diation generation from a system forced by an temporally 
varying external source. It is difficult to solve analytically 
although a solution was given in under the condition 
of Ho = 0 due to the source term independent of time. 
Obviously both the y and z components of the radiation 
have a strength linearly proportional to vq according to 
Eq. ([5]). The two components have a phase displacement 
fixed at 12 which is determined by the one between VyQ 
and Vzo- Therefore, the radiation at uj+ is CP. The rota¬ 
tion of Vo and the radiation will be reversed provided the 
B held sign is changed. These agree with the simulation 
results above. 



FIG. 3. (a) Temporal waveforms of the THz electric fields Ey 
on the axis and (b) the corresponding spectra, where different 
lines correspond to different cjc (uc = ITHz corresponds to 
Bo = 35.7T). 

To further study the radiation features, we vary the B- 
held strength and gas density in the following simulations 
as shown in Figs. 3 and 4. Figure 3 illustrates that with 
an enhanced B held and ojc = 10 THz, the radiation 
frequency, polarization, and waveform agree well with the 
analysis since the condition given in Eq. o is sufficiently 
met. With uJc = 2. THz the simulation results roughly 
agree with the analysis. When uJc = oJp = 1 THz, the 
radiation still has the Ey component and its spectrum 
within 0.5-1.5 THz agrees with Eq. ©■ Its waveform 
attenuates with time, approaching the one without the 
B held, because its frequency is close to Wp. 


















4 



frequency (THz) 

FIG. 4. (a) Temporal waveforms of the THz electric fields Ey 
on the axis and (b) the corresponding spectra, where different 
lines correspond to different ujp or gas densities. The B field 
strength is fixed at 357T. The inset in (a) shows the radiation 
fields Ey and E^ with ojp = 4 THz and the components below 
2 THz are filtered. 


In Fig. 


10^® cm 1.1 X 10^^ cm ^ 


4 we take different gas densities: 4.9 x 
and 1.96 x 10^^ cm“^ with 
the corresponding Wp = 2, 3 and 4 THz. The radiation 
frequencies agree well with Eq. dSD- The radiation has 
a many-cycle waveform but the temporal attenuation of 
the waveform becomes more obvious with the growing 
density since the frequency Cli+ is closer to Wp. The radi¬ 
ation is nearly CP even with Wp = 4 THz, as observed in 
the inset in Fig. 4(a), provided the lower-frequency com¬ 
ponent is filtered. Note that as uj- increases with the 
growing Wp, the field envelope tends to oscillate around 
the axis {Ey = 0). 


TABLE I. Strength (MV/cm) of CP THz radiation of the uj+ 
component as a function of Wc (the row) and utp (the column), 
where ljJc = 2THz corresponds to Bq = 71T. 


strength is nearly not changed with the B field. Multiply¬ 
ing the electron motion equation by the electron velocity 
V, one obtains dv^/dt = —2eE-v/me, where E is the laser 
electric fields. When the B field satisfies Eq. 0 and tens 
of fs laser durations are considered here, the rotation of v 
from the laser polarization plane is slight during the laser 
interaction with the electron and therefore, the net gain 
of the electron energy (also the radiation energy) is nearly 
the same as the case without the B field. With a given 
radiation energy, the radiation strength will decrease lin¬ 
early with its frequency oj+ ~ Wc, i.e., E'^^^ oc l/wc. 
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FIG. 5. [(a) and (c) (Temporal waveforms of the THz electric 
fields on the axis and [(b) and (d)] the corresponding spectra, 
where the gas density is taken as 1.96 x 10^^ cm“® (cjp = 
4THz) and different lines correspond to different Wc (uJc = 
0.5THz corresponds to Bq = 17.8T). 

Next, we consider relatively low B fields with ujc < ojp 
as shown in Fig. 5. EP radiation is also generated even 
with Bq = 17.8T or uic = 0.5THz. The amplitude of Ey 
grows with the B-field strength since more electron en¬ 
ergy is transferred to the y direction by the B field. Both 
Ey and E^ have near single-cycle waveforms since their 


We list the radiation strengths obtained in simulations 
as a function of lOc and ujp in Table I. It is shown that 
the radiation strength roughly follows: 


7?'^ + 
^THz 


OJ" 


OC 


UJn 


( 7 ) 


According to Eq. 0, the strength scales linearly with 
the plasma density or the net current strength, i.e., 
oc Wp. With a given plasma density, the current 


77 '^+ 

^THz 


frequencies separated) and the relative phase is changed 
from O.Stt to around O.Stt as Wc grows from O.STHz to 
4THz. 

The radiation with either a high or low B field is gen¬ 
erated due to the gyrational motion of plasma electrons 
under the B field. This field slightly affects the gas 
ionization responsible for the current formation. Hence 
the magnetic approach can be extended to other laser- 
plasma-based THz emission schemes, which has been ver¬ 
ified by our simulations with an asymmetric-laser scheme 
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(4^ . Besides, we have performed a 3D PIC simulation 
and observed the same result as a 2D simulation with 
the same parameters as in Fig. 1 but with a plane laser 
driver. The result also approaches that with tq = 150/im 
in Fig. 1(c), indicating that our study is valid with a 
larger spot radius tq. 


In summary, we have demonstrated a unique EP or 
CP, narrow-band THz source if a static B field is ap¬ 
plied. With a high B field at a lOOT scale, the radia¬ 
tion shows two frequencies: the lower is nearly dc and 
the higher (central frequency) almost at lOc- Therefore, 
the central frequency can be adjusted linearly by the B- 
field strength. The radiation rotation can also be con¬ 
trolled by the B-field sign. With the B field decreased 
to the lOT scale, EP radiation is still generated but be¬ 
comes broadband and single cycle. To fully apply this 
scheme under the high-B-field condition in Eq. (4) to 
the whole THz band, the B-field strength should be tun¬ 
able within 3.57T to 357T (corresponding to O.lTHz to 
lOTHz). Besides the B-fields generated in traditional 
ways 31-3^, one may em ploy n ovel nanosecond-laser- 
driven superstrong B-fields [3J,1^ with the strength con¬ 
tinuously tunable by the laser energy and with a duration 
of nanoseconds (far longer than THz radiation cycle). 
With this source, one can realize the magnetically con¬ 
trolled two-color scheme in an all-optical way, where the 
B-field generation is easily synchronized with the two- 
color laser driver. 
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